Local-moment formation in gapless Fermi systems 
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Perturbative scaling is applied to the Anderson model for a localized level coupled to a Fermi system 
in which the density of states varies like |e| ? near the Fermi energy (e = 0). This model with r = 1 or 
2 may describe magnetic impurities in unconventional superconductors and certain semiconductors. 
The pseudo-gap is found to suppress mixed valence in favor of local-moment behavior. However, the 
magnitude of the exchange coupling J in the local-moment regime is reduced, thereby decreasing 
the parameter range within which the impurity becomes Kondo-screened at low temperatures. 



75.20.Hr, 72.15.Qm, 74.70.Tx 

There exists a class of "gapless" Fermi systems which 
exhibit a pseudo-gap in the effective density of states 
p(e) at the Fermi level, taken to be e = 0. For instance, 
the valence and conduction bands of certain semiconduc^. 
tors — including Pbi-^Sn^Te at, a critical composition,!!! 
and PbTe-SnTe heterojunctionsrJ — touch in such a way 
that, for small |e|, p(e) is proportional to |e| d_1 in d spa- 
tial dimensions. The quasiparticle density of states in 
an unconventional superconductor can n vary like |e| or 
e| 2 near line or point nodes in the gap.H Heavy-fermion 
and cuprate superconductors are candidates_for this be- 
havior. Electrons in a strong -magnetic fieldu and exotic 
phases of the Hubbard modelH are also predicted to ex- 
hibit a linear pseudo-gap in two dimensions. Finally, the 
single-particle density of states in the one-dimensional 
Luttinger model varies like |e| 2a , where a is a positive 
number which varies continuously with the strength of 
the repulsive Coulomb interactions in the bulk.cl _ ._. 

Recently there has been considerable interestLnl3 in 
the behavior of magnetic impurities in gapless systems 
having a power-law density of states, p(e) — po\e\ r . Poor 
man's scalingTor the spin-^ (impurity c eseneracy N = 2) 
Kondo modelQ and large- N treatmentsoEl indicate that a 
Kondo effect (i.e., complete quenching of a local mag- 
netic moment in the limit of low temperatures T) takes 
place only if the antiferromagnetic electron-impurity ex- 
change poJ exceeds a critical value, pqJ c ~ r; otherwise, 
the impurity decouples from the band. A large- N. study 
of magnetic impurities in gapless superconductorsQ yields 
similar results, except that for r < 1 or N = 2, any fi- 
nite impurity concentration drives J c to zero. Numerical 
renormalization-grojip calculations for the N = 2 case, 
both atEl and nearLH particle-hole symmetry, show that 
J c — > oo for all r > |, while for r < j the strong-coupling 
limit exhibits anomalous properties, including a non-zero 
moment. Away from this symmetry, a finite value of 
J c , roughly proportional to r, is recovered; for J > J c 
the impurity spin is completely screened at T — 0, but 
an electron phase shift of tt suggests that the impurity 
contribution to the resistivity vanishes, instead of taking 
its maximal nassible value as it does in the conventional 
Kondo effect H 

The Kondo model presupposes the existence of a lo- 



cal moment, i.e., an impurity level having an average 
occupancy (rid) = 1. This paper reports the first sys- 
tematic exploration of local-moment formation in gapless 
systems. Poor man's scaling^ is applied to the Anderson 
impurity modeljij in which mixed-valence (0 < (rid) < 1) 
and empty- impurity ((rid) ~ 0) regimes compete with 
local-moment behavior. Concentrating on the case of 
a localized level which lies within a power-law pseudo- 
gap, we show that the reduction in the density of states 
near the Fermi level has three main effects, each of which 
grows more pronounced as r increases: (1) The mixed- 
valence region of parameter space shrinks, and for r > 1 
disappears altogether. (2) The local-moment regime ex- 
pands. (3) The value of the Kondo J on entry to the 
local-moment regime is reduced. Since the threshold J 
for a Kondo effect rises with r (see above), these results 
imply — at least in the cases of greatest interest, r = 1 
and 2 — that there is a large region of phase space in 
which the low-temperature state has an uncompensated 
local moment. This should be contrasted with systems 
having a regular density of states, in which an Anderson 
impurity is always quenched at zero temperature. 

We start with the Anderson modeljij written in one- 
dimensional form: 



H = l de ec\ a c ta + e d n d + U n d ^n d i 

+ J2f de v ^p& ( c - d - + h - c -)- (!) 

„_* , J—D 



<r=T,4' 

The conduction band is taken to be isotropic in momen- 
tum space and to extend in energy over a range ±Z? about 
the Fermi energy; the operators c eo . are normalized such 
that {c\ a , c £V /} = S(e — e')5 CTi<r /. The localized level is 
described by its energy td relative to the Fermi level, and 
the cost in Coulomb energy, U > 0, when it is doubly oc- 
cupied. We have assumed the hybridization between the 
band and the impurity is purely local; V can be taken to 
be a positive real number. The interesting physics of this 
model occurs when V and \ed\ are smaller than U and D. 

Pure power-law density of states. We first consider the 
density of states introduced in Ref. 0: 
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Po \e/D\* 
0, 



\z\<D; 

otherwise. 



(2) 



The exponent r can take any non-negative value, with 
r = representing a constant density of states; pa is 
chosen so that J_ D p(e)de = 1. (Later in the paper we 
will examine a more realistic case, in which the power-law 
variation is restricted to the vicinity of the Fermi level.) 

In order to understand the behavior of Eq. Uu) at low 
temperatures T, we apply poor man's scaling.Ej In this 
approach, electronic states with energies e| 3> T are pro- 
gressively integrated out, yielding an effective descrip- 
tion of the problem in terms of fewer degrees of freedom. 
Consider first an incremental reduction of the bandwidth 
from D to D' = D(l + 5\nD) < D. The aim is to rep- 
resent the same physical system by an effective Hamilto- 
nian of the form of Eq. (|l|) and a density of states given 
by Eq. (||), but with D replaced by D' in both equa- 
tions. This requires the couplings entering Eq. ([!]) to 
be adjusted to account for the states which have been 
eliminated. Provided that 5D, e d , and V remain small 
compared to D, these renormalizations can be computed 
within perturbation theory. The band reduction can then 
be iterated, leading to differential equations for the cou- 
plings as functions of the effective bandwidth. 

Poor man's scaling neglects higher-order corrections 
which conceivably could accumulate to become im- 
portant at low energies. However, the scaling pic- 
ture presented below is supported, both qualitatively 
and quantitatively, bjf non-perturbative renormalization- 
group calculations Jl3 Details of the non-perturbative 
treatment lie beyond the scope of the present paper. 

In previous implementations of poor man's scaling 
for thp-|p^ndcrson model with a constant density of 
stateSjtllHI it was found that eliminating all states with 
energies D' < \e\ < D produces a leading correction 



Se d = -p(D)V 2 5lriD = -(T/w)6\nD, 



(3) 



where T = irpoV 2 . The corrections to U and V enter at 
higher order in small couplings, and can be neglected. 

For a power-law density of states, an additional cor- 
rection is required because the coupling entering Eq. (Q) 
is not V, but y/pJejV = ^T\^~D\ r : . Replacing D by D' 
in Eq. (||) increases p(e) by a factor of (D/D') r . T must 
be reduced by the same factor so that the physical cou- 
pling remains unaffected by the change of variable, i.e., 
we require r | e/D\ r = (T + ST) \ e/D'\ r , which gives 



ST = rTSlnD. 



(4) 



The novel behaviors of an Anderson impurity in gapless 
systems all stem from this correction to T, which has no 
counterpart for a constant density of states. 

Equations (§) and (f|) can be integrated to give the 
effective couplings e d (D) and T(D) at an arbitrary D in 
terms of the initial bandwidth D and the bare couplings 
e° = e d (D ) and T = r(£> ): 
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FIG. 1. Scaling of the impurity energy ed with the 
bandwidth D, shown schematically for (a) r = 0; (b) 
< r < 1. Renormalization of begins on entry into the 
valence-fluctuation (VF) regime, at D w U. Each trajectory 
ends at a crossover to local- moment (LM), empty-impurity 
(EI) or mixed-valence (MV) behavior. Note that the range of 
bare impurity energies that scale to the LM regime (e d < e dt ) 
is much greater in (b) than in (a). 
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r(D) = r 

e d (D)=4+^[ur-iD/D n )'}. 



(5) 
(6) 



[The quantity u = min(l, U/ Do) appears because, strictly 
speaking, Eq. (||) applies only in the range D < U; the 
scaling of is negligible for D > U ] For a constant 
density of states, by contrast, only renormalizestD: 



e d (D) = e° d + (Tq/tt) ln(uD /D) (r = 0). 



(7) 



It should be emphasized that U does not renormalize 
significantly for any r. 

It is instructive to compare the behavior of systems 
with r = and r > as the bandwidth D is progres- 
sively reduced. We assume that the bare parameters are 
such that |e°|,r < U < D . Initially, all four impurity 
configurations-are active, and the impurity susceptibility 
Xi mp satisfies^ Tximp « 1/8. 
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Once D is scaled below U, the doubly occupied im- 
purity state becomes frozen out, and the system enters 
the valence-fluctuation regime. Real charge fluctuations 
between the remaining impurity states lead to a suscep- 
tibility Tximp — 1/6. Here, ed scales with D according to 
Eq. or Eq. ([?]). This renormalization is represented 
schematically in Fig. [j], which shows trajectories ed vs. 
D, starting from different bare values. (The value of To 
is taken to be the same for all curves.) Each trajectory 
in Fig. |l| terminates when either T or \ed\ grows to equal 
D. At this point perturbation theory breaks down, and 
there is a crossover to one of three regimes in which real 
charge fluctuations are frozen out: 

Local moment: If ed becomes sufficiently large and neg- 
ative, the impurity acquires a spin, and the susceptibility 
rises to Tximp ~ 1/4. The scale for this crossover is a 
solution of the equation Dlm = ^^<i{Dlm)] see curves 
(i) in Fig. | 

Empty impurity: If, instead, Cd becomes equal to +D 
[e.g., curves (ii) in Fig. [I]], the localized level is completely 
depopulated and T\imp drops rapidly to zero. 

Mixed valence: Finally, T/D may become of order 
unity, at a bandwidth Dmv = ^{Dmv)] see curves (iii) 
in Fig. |l|. In this case, T\imp — > 0, but even at D = 
the value of (rid) differs significantly from both and 1. 

There are a number of notable differences between the 
cases r — and r > shown in Fig. [l]: 

(1) The trajectories in the valence-fluctuation regime 
are flatter for the system with a pseudo-gap. The larger 
the value of r, the more the renormalization of td is in- 
hibited by the decrease in T. Consider, for example, the 
maximum possible shift in the impurity energy: accord- 
ing to Eq. (||), £rf(0) — e d — (r M r )/(7rr), whereas for 
r = the shift in ed is unbounded. For r > \ (say) 
and r -C |e2|, it is a reasonable first approximation to 
neglect the renormalization of ed altogether. 

(2) As r increases from zero, the crossover scale for the 
mixed- valence regime is pushed down: 
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r • (r /c ) r/(1 - r) , 
o, 



< r < 1; 
r > 1. 



(8) 



For r > 1, the ratio T/D always decreases under scaling, 
which completely rules out mixed-valence behavior; in- 
stead, the system must eventually enter either the local- 
moment regime or the empty-impurity regime. 

(3) The depression of the mixed-valence scale Dmv 
and the flattening of the scaling trajectories both tend 
to widen the range of bare impurity energies which even- 
tually result in local-moment behavior. Curves (iv) in 
Fig. |l| show a case in which the pseudo-gap diverts a tra- 
jectory away from the mixed- valence regime to intersect 
the local-moment line. To quantify this trend, let e d „ be 
the largest (least negative) bare impurity energy which 
flows to the local-moment regime (see Fig. [|). This en- 
ergy is given implicitly by the equation Dmv = Dlm- 
Using Eqs. (^J) and (j^), one obtains 
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FIG. 2. Boundary of the local-moment regime, |e°* |/ro , 
plotted vs. To for different values of r. Any bare impurity 
eventually leads to local-moment behavior. 



energy e° d < e° ds „ 
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(9) 



The plot of e^/To in Fig. || clearly shows the expansion 
of the local- moment region with increasing r. 

(4) At the point where a scaling trajectory enters 
the local-moment regime, the Anderson moidel can be 
mapped via a Schrieffer- Wolff transformation^ onto the 
Kondo model with an effective exchange coupling 
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(10) 



For given impurity parameters (e?, U, and To), the 
Kondo J for r > is reduced compared to that for a 
regular density of states, due both to the depression of T 
and to the weaker renormalization of ed- A lower bound 
on the reduction factor, obtained by neglecting the renor- 
malization of ed, is |e°/£>o| r - This effect is illustrated in 
Fig. ||, which plots poJ as a function of e^, for U = oo, 
r = Do/10, and several values of r. For r < 1, p J 
rises to reach 2 /it (the dashed line in Fig. ||) at e d = e^, 
the boundary of the local-moment region, whereas for 
r > 1, poJ decreases instead. Note that only for r = 
and r = 1/4 is the condition poJ > r for the existence 
of a Kondo effect satisfied over any extended range of e d . 
This observation extends to other values of Tq and U. 

Restricted power-law density of states: In most gapless 
systems, the power-law variation of the density of states 
does not extend over the entire band in the manner as- 
sumed in Eq. (||). We therefore repeat the preceding 
analysis for a more realistic density of states which rolls 
over to a constant beyond a region of width ±A about 
the Fermi level, i.e., p(e) — po\e/A\ r for |e| < A, but 
p{e) = p for A < |e| < D. 
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FIG. 3. Exchange coupling poJ on entry to the lo- 
cal-moment regime, plotted vs. |e^| for different values of r. 
Notice that the exchange coupling decreases as r increases. 



At energies much greater than A, scaling should pro- 
ceed very much as for a constant density of states. It is 
quite possible for the system to pass out of the valence- 
fluctuation regime before the pseudo-gap can have any 
real effect. However, we are more interested in values 
of e° and Tq which are sufficiently small that the band- 
width can be scaled into the range D < A, where Eqs. (§) 
and (|]) must apply. The subsequent renormalization of 
r and £d is identical to that for a system having a pure 
power-law density of states, but with a bare bandwidth 
A, and a bare impurity energy £d(A) calculated from 
Eq. (0). The qualitative effects of the pseudo-gap should 
therefore be the same as those found above, although the 
magnitude of these effects will certainly decrease as the 
width of the pseudo-gap becomes smaller. In particu- 
lar, there will remain an expansion of the local moment 
regime, in which the Kondo J will still be reduced rela- 
tive to the case r — by a factor of at least |ed(A)/A|' r . 

Finally, we note a parallel between our findings for 
an impurity in a non-interacting background and results 
showing that an Anderson impurity in a Luitinger liq- 
uid has an expanded local-moment regime £3 This sim- 
ilarity appears surprising, since in an interacting one- 
dimensional electron gas, independent fermionic ex ci- 
tations are replaced by collective bosonic modes.H It 
seems, though, that bulk interactions affect the valence- 
fluctuation regime of the impurity only by generating 
a power-law spectrum of one-electron states which hy- 
bridize with the localized level. The similarity with 
the non-interacting case does not extend into the local- 
moment regime: there is no threshold value of J-4iequired 
to achieve a Kondo effect in a Luttinger liquid. l3 

In summary, we have investigated local-moment forma- 
tion in Fermi systems having a density of states which 
vanishes as |e| r near the Fermi energy. The pseudo- 
gap strongly suppresses hybridization between conduc- 



tion electrons and an impurity level, and weakens the 
renormalization of the impurity energy. These effects in 
turn expand the range of bare impurity parameters that 
lead to localization of a spin at the impurity site, but 
reduce the value of the Kondo exchange coupling of this 
spin to the band, and thus decrease the likelihood that 
the moment will be Kondo-screened at low temperatures. 
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